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Abstract. We generalize the definition of localization length to disordered systems driven by a time- 
periodic potential using a Floquet-Green function formalism. We study its dependence on the amplitude 
and frequency of the driving field in a one-dimensional tight-binding model with different amounts of 
disorder in the lattice. As compared to the autonomous system, the localization length for the driven 
system can increase or decrease depending on the frequency of the driving. We investigate the dependence 
of the localization length with the particle's energy and prove that it is always periodic. Its maximum is 
not necessarily at the band center as in the non-driven case. We study the adiabatic limit by introducing 
a phenomenological inelastic scattering rate which limits the delocalizing effect of low-frequency fields. 

PACS. 72.15.Rn Localization efi'ects - 73.20. Fz Weak or Anderson localization - 73.21. Hb Quantum wires 
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• 1 Introduction 

■ The seminal paper of Anderson on the absence of difFu- 
] sion in certain random lattices started a new era in 

■ the study of the effects of disorder in condensed matter 
physics. Real materials always contain a certain degree of 

, disorder, since the atomic structure is never perfectly reg- 
' ular, and many physical properties are either influenced 
I or even mostly determined by this randomness 2 . As the 

■ translational symmetry is broken, electrons in disordered 
' lattices are not delocalized Bloch waves and can be spa- 
, tially confined. One of the main quantities of interest in 

disordered systems is the localization length A of the elec- 
' tron wave- functions. A system of typical size L > X will 

■ behave as an insulator while a system with L < X will 
I behave as a conductor |2I3| . 

The Anderson model in one dimension has been used 
as a fundamental tool for the understanding of the prop- 
, erties of quantum and molecular wires since the beginning 
1 of the field of molecular electronics @] . In one dimension 
it can be shown in a mathematically rigorous way that 
all eigenstates of the Anderson model are exponentially 
localized in the thermodynamic limit 

On the other hand, the use of a time-periodic driv- 
ing field in nano-technological applications has emerged 
in recent years as a way to control the properties of possi- 
ble future devices |H] and as a source of many interesting 
new effects [J]. For example, in semiconductor superlat- 
tices in the presence of THz radiation, absolute negative 
conductance and dynamical localization have been ob- 
served experimentally j8l9| . In the field of molecular elec- 
tronics, a ratchet-like effect was shown experimentally in 
the photocurrent in a self-assembled monolayer of asym- 



metric molecules IIUI. In terms of future applications, it is 
very important to understand the effect of periodic driv- 
ing in the properties of quantum wires. It can be argued 
that the localization behavior of this kind of systems is its 
most basic property. It is the goal of this paper to give the 
basic foundations for the extension of localization theory 
to periodically driven systems and to study the prominent 
features of a simple model that includes both disorder and 
periodic driving. 

In the absence of disorder, and for a 1-D tight bind- 
ing system with a band width A, it was shown that, for 
high-frequencies, the driven system is equivalent to a non- 
driven system with a renormalized band-width |15ll'5| , 



A 



A Jo 



(1) 



where Vac is the amplitude of the ac-driving field and d 
is the spatial period of the lattice. At the zeros of this 
Bessel function the width of the band goes to zero, the 
group velocity of an electron wave-packet becomes zero 
and the electron becomes effectively localized (however, 
the Floquet eigenstates continue to be extended). This 
strong localization of a particle due to the effects of a 
periodic driving has been called dynamical localizationfT^ 
or coherent destruction of tunneling '14'. 

The effect of a harmonic driving in the localization of 
a disordered one-dimensional system was first studied by 



Holthaus et al. |llll2j , who showed that a driving field can 
have strong effects in the localization of a particle moving 
in a disordered potential. When disorder is introduced it 
is well-known that all the states in the one-dimensional 
lattice are localized, with their degree of localization con- 
trolled by the ratio between the bandwidth A and the 
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strength of the disorder W. Holthaus et al. have argued 
that in the high-frequency hmit the Floquet states them- 
selves become locahzed (in contrast with dynamical local- 
ization) and their localization should depend on the ratio 

In a recent work | l/Zj we were able to generalize the 
definition of localization length using a Floquet-Green 
function formalism. Through a numerical calculation of 
the Floquet-Green function for a one-dimensional tight- 
binding model with diagonal disorder plus an external 
dipolar driving field we could confirm quantitatively that 
in the high-frequency regime the effect of the driving could 
be seen as a renormalization of the band-width. More 
interestingly, in that work it was also shown that low- 
frequency driving can reduce the localization of the elec- 
tronic Floquet-states. 

It is the purpose of this paper to present, in a de- 
tailed and self-contained manner, the theory first used in 
|17| . We also present new analytical calculations for the 
Floquet-Green function in the high-frequency regime and 
explore additional features of our model regarding the pe- 
riodic dependence of the localization length on energy as 
well as the effect of inelastic processes in the delocalizing 
effect of low- frequency driving fields. 

Section[21is devoted to the Floquet-Green operator for- 
malism which is the basis of all our results. Here we show 
its relationship to the transport properties of driven sys- 
tems. The generalization of the definition of localization 
length is discussed in section |21 In section 21 we introduce 
the Hamiltonian and the method that we follow to ob- 
tain the localization length for this system. In section O 
we examine the behavior of the localization length for dif- 
ferent amplitudes of the driving field and in different fre- 
quency regimes. In particular, the results for A in the case 
of high-frequency driving will be studied in subsection l5.ll 
In this section we also derive some new analytical results 
for driven non-disordered systems which explain some fea- 
tures of the localization length of a disordered system in 
the high-frequency limit. Some previous obtained results 
for low-frequency delocalization will be presented in sub- 
section New results for A as a function of energy will 
be shown in sectional The previously unexplored effect of 
an effective inelastic scattering rate of the electrons on the 
delocalization properties of low-frequency driving will be 
studied in subsection [71 Finally, in section |S1 we give some 
concluding remarks and perspectives. 



2 Floquet-Green operator for time-periodic 
systems 

The pioneering work of Shirley JHIj Zel'dovich ul9_ and 
Sambe laid down the theoretical foundations for a 
complete treatment of time-periodic potentials, based on 
the same mathematical tools already developed for time- 
independent potentials. Of significant importance among 
these tools is the Green function. A Floquet-Green func- 
tion method for the solution of radiative electron scatter- 
ing in a strong laser field was introduced by FaisaH24l. 



The definition and application of the Green function for- 
malism that fully exploits the periodic time-dependence 
of the Hamiltonian has not been done until recently. In 
this section we show the complete Floquet-Green opera- 
tor formalism for general time-periodic Hamiltonians, in 
the way that was introduced by one of the authors in pre- 
vious works 22 23 . 

We start by considering a general Hamiltonian of the 
form: 

H{t) = iJo + 2Vcos{ujt), (2) 

were Hq and V are Hermitian operators in the Hilbert 
space {TL) of the system. Because of the periodicity of the 
Hamiltonian, according to Floquet 's theorem, the solu- 
tions to Schrodinger's equation ih-^\^{t)) = H{x, t)\^{t)) 
are of the form 

\<lr-{t)) = e--*/'^I(/.^(i)), (3) 

where \r{t)) ^ irit + ^)). ^ 

Inserting this into Schrodinger's equation one arrives 
at the eigenvalue equation 

H^itWit))^e\rit)), (4) 
where H^{t) is defined as 

H^{t)^H{t)-ih^^. (5) 

As pointed out by Sambe [201 j since Eq. I^J is an eigen- 
value equation, it can be solved using the standard tech- 
niques developed for time- independent Hamiltonians, pro- 
vided we extend the Hilbert space to include the space of 
time-periodic functions. In this extended space, the time 
parameter can be treated as another degree of freedom 
of the system. A similar concept is used in classical me- 
chanics, and gives rise to the concept of a " half degree of 
freedom when dealing with time-dependent Hamiltonians. 

A suitable basis for this extended Hilbert space {TZ) is 
{|a} ® |n),...}, where {|a),..} is a basis for the Hilbert 
space Ti. of the system, and we define {t\n) — e"™"^*, 
with n integer. Clearly {|n),..} spans the vector space 
(T of periodic functions, and therefore, TZ = Ti^T. In 
this basis, Eq.Q becomes a matrix eigenvalue equation 
of infinite dimension with an infinite number of eigen- 
values. It is not difficult to prove that if Ca.p is an 
eigenvalue with corresponding eigenvector (t)), then 
ea^p + mhio is also an eigenvalue with corresponding eigen- 
vector |0eo.P+™'i'^(t))= e""'^*|0^°.p(i)}. Accordingly, the 
eigenstate corresponding to the eigenvalue Cq + mhuj has 
the same structure as the eigenstate corresponding to ea,p, 
except that it is displaced by mhuj on the energy axis. 
Because of this, to find all the eigenvectors and eigenval- 
ues of the Floquet Hamiltonian one needs only to con- 
sider — < e < ^huj. We will use the letter e to refer 
to the Floquet eigenvalues restricted to this interval and 
call them " quasi-energies" . Clearly, any Floquet eigen- 
value Ba^p can be written as Ca.p — Sa + pfu^ for some 
— •^hjj < Ea < \^ a-nd some integer p. It can be shown 
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that, in general, there are N distinct quasi-energies (ex- 
cept for accidental degeneracies) if the Hilbert space H is 
N-dimensional. 

This periodic structure in the eigenvalues does not 
mean that the "sideband" eigenstates have no relevance; 
they are also valid solutions of Eq. and are essential 
for completeness in the extended Hilbert space TZ 21 . 

The Floquet-Green operator corresponding to Eq. 
is defined by the equation (see |22I23| '). 



[IE - H^{t')]G{E,t', r) = I5r{t' - t") 



(6) 



where 5r{t) is the r-periodic delta function (t = ^) and 
I is the identity operator in H, which we will omit from 
now on. 

Notice that this equation differs from the usual time- 
dependent Green function on two instances: The Hamil- 
tonian H{t) is replaced by H^{t), as defined in Eq.©, 
and also S{t) is replaced by Sr{t). With this definition, the 
properties derived from the periodicity of the Hamiltonian 
and its eigenfunctions are built-in features of the Floquet- 
Green function. In terms of the complete (infinite) set 
{|(/)"'P(t))} of eigenfunctions of the Floquet-Hamiltonian 
fEas l4l5|l . the solution for Eq.® can be written as 



(7) 



^From the previous discussion about the eigenvalues 
and eigenfunctions of the Floquet Hamiltonian, we can 
write the Floquet-Green operator entirely in terms of the 



eigenfunctions |^"'°(i)), which correspond to values 



between —\fy-^ and \huj : 



G{E,t\t") = Y,Y.^ 



E — Ca — 



(8) 



where 7 — l,...N for Ti being N-dimensional, and p 
—00, 00. 

Operating on both sides of this equation with 

J_ (-T 
r2 Jo Jo 



1 r r g*ma;t'g-»n^t"^^./^^. obtain 



G"'"(£;) = J2 



1 



\<l>nf+n){K+v\: (9) 



where 



G""'"(£;) = ^ 



Jo 



.^rnu^t' ^-^nut" G{E,t' ,t")dt" dt' 



and 



1 
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At this point, a technical comment is in order. Notice 
that the parameter E in the Floquet-Green function is not 
a quasienergy since it is not restricted to the first Bril- 
louin zone. The Floquet-Green function has poles along 
the whole real axis. When considering a closed system 
with a Hamiltonian periodic in time, it is well-know that 



the energy is not a conserved quantity and that the sys- 
tem state is specified by its quasienergy. From this point 
of view, one might think that E should be restricted to the 
first Brillouin zone. However, consider the physical situa- 
tion in which the particle comes from outside the system 
and one is interested in the probability for the different 
processes that can occur to this particle in the course of 
its interaction with a potential. These quantities can all 
be obtained from the Green fmiction, and in the case of 
time-periodicity, from the Floquet-Green function. In par- 
ticular, G'™^"(£') gives the probability amplitude for a pro- 
cess where the incident particle with energy E enters in 
(Floquet) channel n, and leaves the system through chan- 
nel m. The meaning of these "channels" is clear from the 
following property of the Floquet-Green function: From 
Eq.lUni), 



G"'"(£;) = G' 



i — k.n—k 



{E + khw). 



This means that the energy of the incoming particle can 
be thought of as determining the incoming and outgoing 
channels. In fact, G™""{E) can be interpreted as giving the 
probability of a particle that comes in with energy E + n, 
absorbs m — n photons and leaves with energy E+m. From 
this it is clear that to describe an incident particle with 
energy E one can either decompose E = e + nhuj and use 
the quasi-energy e along with an incident channel n, or 
simply use E as the energy and assume that the incident 
channel is n = 0. In the following, we adopt the last view 
which we believe is the most natural one. 

^From the previous discussion, one can see that the 
quantities G^'^{E) provide all the information of the 
driven system. 



G'''°{E) = 



1 



E 



(11) 



It is easy to show that in the hmit V ^ 0, G°'^{E) 
G{E), where G{E) is the usual Green function for the 
autonomous system. Also, in this limit, G'^'^°'°(i?) 
(see rcf. 23^ for an equation that gives G^'^{E) in terms 
of G"'*' (£')). The Floquet-Green operator components 
G^'^{E) are important because transport properties of 
driven systems have been formulated in terms of these 
components, which play a role similar to the Green oper- 
ator in the Landauer formalism for conduction More 
specifically, it was found that the average current through 
a 1-D driven system (coherent regime) can be expressed 
as 



/ = 



OO „ 

jdE {t^^\e)UE) - T^',\E)f,{E)} , 



(10) where 



A:— — OO 



T^^'{E) - FeiE + khu;)rr{E)\Gl'l{E)\ 



fe,0/ 



(fe)/ 



.^,>{E)=rr{E + khLo)r,{E)\G):°,{E)\\ (12) 

denote the transmission probabilities for electrons from 
the right(left) lead respectively, with energy E and fi- 
nal energy E + khu, i.e. the probability of a scattering 
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event under the absorption(emission) of \k\ photons, if 
k > 0{k < 0). This expression for the average current 
converges to the well-known Landauer-Biittiker formalism 
|26| in the limit when the driving amplitude goes to zero. 



3 Definition of localization length for driven 
disordered systems 

The localization length for 1-D disordered non-driven sys- 
tems has been defined in terms of Green functions [H], 
through the well-known relation 



HE) 



hm -(ln|GiL(£;)|). 

^— >00 J_j 



(13) 



This definition makes use of the proven fact that in a dis- 
ordered 1-D potential, the wave functions for any energy 
decay asymptotically in an exponential way with distance: 

Using Ea. (|13|l . and for a tight-binding Anderson 
Hamiltonian (band-width A and diagonal disorder on 
the on-site energies uniformly distributed in the interval 
[—W/2, W/2]), the localization length is known to behave 
as 127] 

W. 



A = 6.56 



(14) 



valid for A/W not too large, and for an energy in the 
middle of the band (£' = 0). Also, it is known that in au- 
tonomous systems, A as a function of the energy E follows 
an inverse parabolic law with a maximum at i? = 0. This 
is so unless \E\ > (Z\ + W)/2 where the density of states 
is null. 

Clearly, in the presence of a periodic driving, the def- 
inition in Eq. (|13() must be modified. Since the Floquet- 
Green operator is the natural extension for time-periodic 
systems of the Green operator for autonomous systems, 
and in view of the expressions for the transport proper- 
ties of driven ID-systems, Ea. (|12|l . we consider the fol- 
lowing quantities as possible extensions of the concept of 
localization length for driven systems: 



\'-{E) L^loL 



In 



G\t{E) 



(15) 



in close analogy to Ea. (|13|l . The different quantities 

G'lf{E) are associated with the probability of a process 
where an electron starts with an energy E at site 1 and 
ends at site L with energy E + kfku. In principle, there 
could be a different localization length associated with 
each one of these processes. However, since the asymp- 
totic behavior of the Floquet eigenstates is exponentially 
decreasing (as we will show next), one can see that in the 
limit L — > cxD, the dominant term in the sum over p in Eq. 
Hll(l always appears {p can be positive or negative) for any 
value of k. This implies that all the quantities G{^{E) de- 
cay at the same rate with L (even though their values can 
be very different, depending on V and uj.) Consequently, 
the quantities A*-*^^ [E) are all identical for V ^ 0. However, 



as mentioned before, in the Hmit F — > 0, G^'°{E) G{E) 
and GC^'^O'") {E) 0, where G{E) is the Green function of 
the time-independent system. This implies that the best 
choice for defining the localization length of a driven sys- 
tem is A'^°^(£'), since it converges to the localization length 
of the autonomous system when V ^ 0. Due to these con- 
siderations, we finally define the localization length for a 
driven system as: 



lim 



In 



G'itiE) 



(16) 



This definition of localization length seems the most 
straight forward and natural way to generalize this con- 
cept for a driven system. It is the connection between our 
localization length and the transport properties of a driven 
system that makes our definition interesting and useful, 
with clear applications and consequences in possible ex- 
perimental setups. 



4 Model and method 

In order to investigate the effect of a driven potential in 
the localization properties of one-dimensional disordered 
systems we have chosen the Anderson Hamiltonian with 
diagonal disorder. The on-site energies £j are distributed 
uniformly in the interval [—W/2, W/2], where W measures 
the strength of the disorder. The driving potential is due 
to the presence of a time-periodic spatially uniform field, 
e.g. the interaction between an electron and an EM wave 
(dipolar approximation) incident perpendicularly to the 
lattice, with its electric field polarized along the lattice 
direction. (From now on we will use a system of units in 
which h = 1, therefore, to will refer to the driving fre- 
quency measured in energy units.) Accordingly, 



H= E,(U + i)01 + lj)01) + 



(17) 



In this work we will always use A = 4. For the calcula- 
tion of G'^''^{E) we use a method developed by one of the 
authors (see for details). For a periodic Hamiltonian 
of the form H{t) = Ho + 2V cos^ivt), where Hq and V are 
any time-independent operators in the Hilbert space of the 
system, the Floquet-Green operator components satisfy 



{E + ku}- Hq)G^'^ - V{G 



ife+1,0 



G 



Sk.o- (18) 



These equations can be solved using matrix continued 
fractions. For the case fc = 0, one gets 



G"^%E) = {E-Ho-V,ffm-\ 



(19) 



where 



with 



V,ff^V+f{E) + V-f^iE), 



V.lfiE) = V- 



-V 



E±luj-Ho-V- 



E±2uj- Ho -V V 



-V 



(20) 
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Fig. 1. (in G''-|'^(0)) as a function of the length of the sys- 
tem L for different values of disorder W , frequency oj, and 
field amplitude v. High-frequency results, uj = 20, show regu- 
lar oscillations on top of the exponential decay, according to 
the equation G''j'^(0) = AJo{2vL /uj) ey.^~'"^^ . Low frequency 
curves do not show these oscillations. The value of A is the 
negative inverse slope of the linear fits. The ensemble average 
was made with 1000 different realizations of disorder. 



The convergence of equation (|2()|l is system specific. 
For our Hamiltonian, Eq.(l), the number of bands neces- 
sary to ensure convergence increases linearly with vL/uo. 
The numerical performance of our method is determined 
by the speed in the calculation of an L x L-matrix inverse 
for each Floquet sideband. 

To obtain the localization length from the Floquet- 
Green function, we plot the ensemble average oflnG°f (0) 
as a function of the length of the system L. In FiglJwe 
show some examples of the results obtained. As expected, 
a straight line fits the data very well. The wave-functions 
decay exponentially with distance. The negative slope of 
these curves corresponds to the inverse of the localization 
length when L :$> X. 



5 Frequency and amplitude dependence of 
the localization length 

5.1 High-frequency localization 

We now show results for the high frequency regime. Some 
results for X(E =^ 0) in this regime were already reported 
in jl7| . The high frequency limit can be characterized as 
the regime in which the absorption or emission of any 
number of photons would leave the particle with an energy 
outside of the region where the eigenenergies of the non- 
driven system concentrate. This region is well known to 
have a width A + W, and therefore, at the center of the 
band this condition is satisfied when ui > {A + W)/2. 

We will first look at the high-frequency regime without 
disorder. For this case one can derive an analytical expres- 



sion for the Green function of the driven system, following 
the work of M. Hohhaus and D.W. Hone [d. We wiU use 
this expression to infer, for this frequency regime, the form 
of the localization length for the disordered case. 

The Floquet eigenstates can be written in the form of 
Houston states pS] . 

I0"'°(i)> = E l^)e-^P (^^1kit)ld - I dT[E{qk{T)) - e{n)] 

(21) 

where qk{t) = k + 2vsaiujt. Here k is the quasimomentum 
of the Floquet state, which is connected to its quasienergy 
through the dispersion relation 

e{K) ^ Jq{2v/ljj)E{h) mod uj , (22) 

and E{k) = — (Z\/2)cos(k). The quantity E{qk{T)) can be 
written in terms of Bessel functions as 

E{qk{T)) = -^Y. -Wcosin + roJT). (23) 

r 

^From this analytical expressions for the Floquet 
eigenstates one can construct the Floquet-Green opera- 
tor, using Ea. pi|l for fc = 0. For that purpose we need 
to calculate the fourier components of the Floquet eigen- 
states, 

|</.;^°) = i rexp(*mc.t')l0"'°(OK • (24) 
Jo 

^From Eq.|2lJ we get, 

1 

i;)- / dt'e^p {imut' + iqk(t')ld 

t' (25) 

~i [ dt"[E{qk{t"))~e{K)]}. 
Jo 

At this point we can use the fact that in the high fre- 
quency regime, A lo and therefore in this integral, to 
lowest order in A/uj, we can ignore the contribution from 
the terms that contain E{qk{T)) and £{hi) (both are pro- 
portional to A). From this the above expression reduces 
to 



|0«;O) = i VlOe""" / dt'e4"'^*'+i5^''™M')] , (26) 

^ ; "'0 

and using the identity 

exp [izsincj)) — exp {ikcj))Jk{z) , (27) 

r 

one arrives at 

IC")=E^"'I0"^-'"(2^V^) ■ (28) 
I 

This expression is only valid at frequencies much higher 
than the band- width A. From this, and using Eq.l^TJ, 
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Fig. 2. High frequency results. Results for different values of 
disorder W and uj. The corresponding results for a non-driven 
system with renormalized bandwidth obtained from numerical 
data are shown with full lines {W = 5) and dashed lines W — 
10. 



the final expression for the quantity G^'j^ for the system 
without disorder, in the high frequency regime is 

1 g-«K('-l) 

^29) 

As can be seen in Fig. ^ for the second and last 
sets of data (w = 20) there are dips in the function 

^ln(|Gj^(0)|)y The thin (red online) continuous line rep- 
resents the function ln(AJo(2ti/ti;) exp(— L/A)), which fits 
the high-frequency data very well. The dips are clearly due 
to the zeros of the Bessel function. These Bessel function 
factors, which appear in the Floquet-Green function of the 
system without disorder, remain there after the disorder 
is introduced and the ensemble averaging is taken. This 
is due to the fact that they do not depend on any of the 
parameters of Hq . 

^From Eq.(j2n|) above, we can also see that, apart from 
the Bessel function factors, the net effect of the high- 
frequency driving is to renormalize the band- width of the 
non-driven system, i.e. instead of -E — E{k) in the de- 
nominator, one gets E — £(«) = E — Jo{2v/uj)E{k). This 
renormalization in the bandwidth was shown by Holthaus 
and Hone to carry into the localization length of a de- 
fect placed in the driven lattice This means that 
A depends on the ratio between the defect energy and 
the renormalized bandwidth. ^From this, one can expect 
that, for a driven disordered system in the high-frequency 
regime, the results for A can be obtained through 

X{v,A) ^ X{0,AJo{2v/uj)). (30) 

As Fig. |2] shows, for different values of uj and for disor- 
der W = 5 and W = 10, the numerical data is in excellent 
agreement with this hypothesis. The minima of A corre- 
spond to the zeros of the Bessel function Jq{2v/uj). The 



slight deviations are more significant for smaller values of 
w, possibly due to next order corrections in v/uj, since we 
expect Ea. (|30|l to be strictly valid only for infinite fre- 
quency. 



5.2 Low-frequency delocalization 

As it was already found in 17 , the behavior of A for low- 
frequency driving is very different from the high-frequency 
case. As it can be seen in FigEl a), A as a function of 
v/uj initially increases. This delocalization can only occur 
for a; < (Z\ -I- W)/2. As a function of w, the localization 
length in the low-frequency regime increases, then reaches 
a maximum and finally decreases. This maximum value 
increases as the frequency is decreased. 

In 17^ an intuitive interpretation of this result was 
proposed. The driving allows the electron to exchange en- 
ergy with the external field and to explore new regions of 
the phase space. For an electron of initial energy £'o, turn- 
ing on an ac-field allows it to exchange energy with the 
field (quantum-mechanically) through an integer number 
of photons. Due to this, the electronic state becomes a 
superposition with energies En = Eq ± nui. Each electron 
path corresponding to each one of these energies behaves 
as an effective channel with a relative weight that depends 
on the value of the field amplitude v. If the first of these 
new channels is outside the band-width it will have a neg- 
ligible probability of having a localization length higher 
than the original one without driving. However, if this 
new channel is inside the band-width of the disordered 
wire it will have a finite probability of having a higher 
localization length than the original path of the electron 
with energy Eq. As the total localization length of the 
system will depend on the less localized of the electron 
paths, the ensemble-averaged localization length will in- 
crease with respect to the autonomous case. Given that 
the number of effective channels that fit into an energy 
region of size Z\ -I- is inversely proportional to the fre- 
quency, the localization length will therefore increase with 
decreasing frequency as it was found in [T7j . 

These findings can have an important impact on the 
transport properties of disordered nanowires in the pres- 
ence of electric ac-fields. In principle, at high-frequency, 
and due to dynamical localization, the ac-field will de- 
crease the electric current through the device, whereas, at 
low-frequency, an ac-field can be used to increase the local- 
ization length so that it becomes greater than the length 
of the wire, thereby producing an increase of several or- 
ders of magnitude in the current that passes through the 
nanowire. 



6 Energy dependence of the localization 
length 

As we have mentioned previously, A as a function of 
the energy for the autonomous system follows an inverse 
parabolic law with a maximum at the band-center E = Q. 
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Fig. 3. Localization length A as a function of the energy for 
W = 5, and CO = 20. 



The aim of this section is to study how this dependence 
is modified in the presence of driving. First of all we will 
show that due to the structure of the Floquet-Green func- 
tion the localization length X{E) is a periodic function of 
the energy with period lu, as it can be seen in Fig.|(31l. It is 
easy to show rigourously that for any V ^ 0, the localiza- 
tion length is periodic on the energy. Before proving this, 
we emphasize that this is not a trivial statement. This 
periodicity is not present in the Floquet-Green operator 
from which the localization length is obtained. This can 
be shown as follows: 

^IVom Ea. H19|) we have 

G'^'^iE + Lj) = {E-Ho-V+f{E + uj)-V-f^{E + u))-\ 

(31) 

and from Eg. lf^ . 

V+f{E + uj)^E + u;-Ho-V [V+^iE)]-' V, (32) 

V-fAE + uj) = V i T^V. (33) 

E-Ho-V-j]{E) 

Using this in Ea. l|31(l . one gets 

GO-"(ii; + c.) = y-iv;+/i?)G"-°(i?)i/+/i?)i/-i 

^ 1 (34) 

+ V-^V+f{E)V'\ 

^From this equation it is clear that in general G^'^{E + 
uj) ^ GO'O(S). 

We now proceed to show that, despite of the previous 
considerations, the localization length is indeed periodic, 
\{E) = X{E + uj). As we have commented before, from 
Eq.lini) it is easy to see that G"'-''{E) = G™~"^°(.E-hnw), 
and in particular, 

G'''~\E + uj) = G^^'^{E). (35) 

To continue our proof, we consider the quantities \^'^ {E) 
defined by Ea. p5|l . but with the upper index in the 



Floquet-Green function replaced by j. From Eq.JHl we can 
see that exchanging the two upper indexes in X^'^{E) we 
get the same value. In particular, \^'°{E) = \°'^{E). Also, 
as we have shown before, all these different localization 
lengths are equal to each other. Specifically, 

A"'°(ii;) - Ai^O(ii;) , 

and 

X'^^'^{E + uj) = \'^'-\E + uj). (36) 
Using this, along with Ea. l35() we arrive at 

\°'^{E) = X^'^E) = X"--\E + uj) = A"'"(£; + uj), (37) 

as claimed. 

In Fig. 13 we can see the clearly periodic behavior of A 
as a function of the energy, for the high-frequency regime. 
In this case, the form of the function X{E) is similar to 
the form for the autonomous system. The only differences 
are the previously explained renormalization of the band 
width and the fact that this structure is repeated period- 
ically. In the low-frequency case the form of this function 
changes. There is not always a maximum at the band cen- 
ter. In this case one can often find a minimum at -E = 
and two maxima located symmetrically around this point 
(due to the symmetrical distribution of disorder used), 
which are separated by a difference in energy equal to the 
frequency of the driving field. In Fig. 0] we can see how 
increasing the amplitude v we can go from the case of the 
two symmetric maxima ai E = —lj/2 and E = u) /2 to the 
band-center maximum at E' = in a case with disorder 
W = b and frequency lu — 2. There is one intermediate 
critical value of v where A(£') is completely flat, v = 2.59 
in the figure. This value coincides with the first local min- 
ima of A as a function of u. In general we can say that 
in the very high-frequency regime there is always a maxi- 
mum for A at = and that for very low-frequency the 
localization length as a function of energy changes very 
little, with small oscillations around an average value. For 
values of the frequency in between these two extremes, the 
behavior of A with energy depends on the specific value of 
the different parameters of the system. 



7 Effect of an inelastic scattering rate and 
the adiabatic limit 

In real systems, inelastic effects will occur in the system 
and the exchange of energy between the electrons and the 
periodic external field will be hampered. In this section we 
want to address the effect of an inelastic scattering rate 
which will appear in physical systems due to scattering 
with phonons in the crystal or by other mechanisms. The 
microscopic modelling of these inelastic effects is out of 
the scope of this work. We introduce here a phenomeno- 
logical parameter, the inelastic scattering time Ti„. We are 
interested in studying the effects of this inelastic scatter- 
ing time in the localization properties of disordered one- 
dimensional systems. We will show next that this inelastic 
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Fig. 4. Localization length A as a function of the energy for 
W = b, frequency w = 2 and three different values of the field 
amplitude v. 



time creates a frequency cut-off for the delocalization in- 
duced by low-frequency driving. We will also address the 
relationship between these results and the adiabatic limit 
w ^ 0. 

We introduce now a phenomenological inelastic scat- 
tering rate F = l/ri„ just as an imaginary part of the 
energy [^Hl in the Floquet-Green function, Eqs. H19|l and 
(|^ . We have studied four cases for disorder = 10. The 
localization length as a function of v/oj is shown in Fig. [S] 
for a)r = 0, b) r = 0.01, c) T = 0.1 and d) T = 0.5. 

The value of A for w = is also reduced due to the 
inelastic scattering rate from A — 1.41 when F — 0.0 
to A = 1.37 when F = 0.1, to A = 1.19 when F = 0.5. 
For F = 0.01 the reduction could not be seen with the 
ensemble average performed. In the high-frequency case 
the behavior of A as a function of v is unchanged by the 
presence of inelastic scattering, apart from this global re- 
duction. We can clearly see the effect of dynamical local- 
ization previously discussed and the shape of the Bessel 
functions when w is bigger than the band-width. 

When we enter the low-frequency regime, the effect 
of the inelastic scattering is to reduce the delocalization 
induced by the low-frequency driving. The effect is com- 
pletely suppressed for frequencies ui < F. See, for exam- 
ple, the case uj = 0.1 when F = 0.5. This can be easily 
understood by considering that since the exchange of pho- 
tons with the external field allows the electron to explore 
more parts of the phase space and to "choose" the best 
ones for propagation, then clearly, energy dissipation due 
to inelastic processes will limit this exploration therefore 
degrading the ability of the driving field to delocalize the 
electron wave function. 



The inelastic scattering rate F is also very important 
to correctly obtain the static limit a; = from our for- 
malism. In the adiabatic limit, when cj — > 0, it has been 
shown by Moskalets and Biittiker 30 , that the Floquet 







Fig. 5. Effect of an inelastic scattering rate in the localization 
length. A as a function of v/ui for disorder = 10 and different 
values of the frequency lj. In a) _r = 0, b) _r = 0.01, c) F = 0.1, 
d)r = 0.5. 
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Fig. 7. Behavior of Xmax as a function of l/uj for W = IQ 
and different values of the inelastic scattering rate F. 



Fig. 6. Adiabatic limit for ln|GiL| as a function of v for a 
particular realization of disorder with W = 10 and F = 0.5. A 
comparison is made with the results for lu — 0.01. 

scattering matrix agrees with the time average (over one 
cycle) of the stationary scattering matrix. This property 
implies that the Floquet-Green function G^{E + iF) in 
such limit will converge to the time-average of the Green 
function of the stationary system. 

1 F 

G''°\E+ir) = - / Go{v COS Lut,E+ir)dt when lo ^ 0, 
T Jo 

(38) 

where Go{v,E) is the Green function of the Anderson 
model with a linear potential. 

We need to include a small imaginary part to the en- 
ergy to assure the convergence of the time integral, since 
without it, the poles of our Green function lie directly 
on the real energy axis. The adiabatic regime, where the 
right-hand side of the above equation is a good approxi- 
mation of the left-hand side, is achieved when uj < F. This 
energy relaxation therefore becomes a limiting time-scale 
for the delocalization effect of low- frequency driving. From 
the numerical data in Fig. Q , we see that the monotonic 
increase of A with 1/u) stops at to > 20F. When the con- 
dition UJ < F is achieved, the value of A is obtained from 
the time-average of the Green function of the static sys- 
tem, which implies that in this regime there can not be 
any delocalizing effect due to the field since in the static 
limit the electric field does not increase the localization 
length for any field amplitude {Xmax = Aq). To obtain the 
correct static limit using our formalism and for the purely 
theoretical case when there is no inelastic scattering, one 
needs to take the limits — > 0, w ^ 0, provided lo < F. 



8 Conclusions 

We have explored non-perturbatively, both for high 
and low frequency, the localization of disordered one- 
dimensional tight-binding lattices using a Floquet-Green 
function formalism which makes use of matrix continued 



fractions 1221 ■ This formalism allowed us to calculate di- 
rectly the localization length (A) generalizing the stan- 
dard definition for autonomous systems. This quantity is 
very important in the study of disordered systems and 
its properties are well known for disordered non-driven 
systems. In a 1-D tight-binding model with diagonal dis- 
order we have found that in the high-frequency limit, A 
is renormalized as a consequence of dynamical localiza- 
tion. In this regime, the localization length always de- 
creases with the amplitude of the driving field. Our re- 
sults show that the high-frequency regime in this model 
is reached for uo > [A + W)/2. For low-frequency, where 
u; < (Z\ -|- W)/2, we found the opposite behavior: A can 
be significatively increased in the presence of a driving 
field. According to this, one can say that Anderson local- 
ization in this kind of disordered system is enhanced by 
high-frequency driving and diminished by low-frequency 
driving. For low frequency driving, each additional Flo- 
quet channel created by the driving provides the electron 
with additional paths. These will differ in their degree of 
localization, with some having a smaller and others hav- 
ing a greater localization length as compared to the non- 
driven case. Since this quantity is given by the paths which 
have maximum extension in space, we conclude that af- 
ter the ensemble average has been performed, additional 
propagation channels should always contribute to increase 
the localization length. This simple picture however, does 
not apply to the high-frequency case because new paths 
introduced by the absorbtion or emission of one or more 
photons always have localization lengths smaller than in 
the non-driven case. 

We have shown that due to the Floquet structure of the 
states, the localization length in a driven system is peri- 
odic in energy, with the periodicity given by the frequency 
of the driving field. For frequencies well above {A + W)/2, 
the maximum of X{E) is always at i? = as it is in the 
autonomous system. For frequencies well below this value, 
the localization length does not change very much with 
energy. For intermediate values of frequency, the localiza- 
tion length can have either a maximum or a minimum 
a.t E = 0, depending on the specific parameters of the 
system. 
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We have explored the limit of very-low frequencies. We 
have shown that the adiabatic limit for this theory is only 
well defined if we add an imaginary part to the energy in 
the Floquet-Green function. This quantity has the physi- 
cal meaning of an energy dissipation rate or inelastic scat- 
tering rate F and has important physical consequences. 
The dynamical delocalization effect discussed in our work 
will be limited by F and will be completely suppressed for 
frequencies to < F. 

The different results analyzed in our work should have 
important experimental consequences in the field of coher- 
ent transport. They lead to new avenues for the control of 
the localization and transport properties of quantum wires 
and also for atoms in optical traps, which have also been 
proposed as a testing ground for Anderson localization. 
Recently, several experimental results with BEG atoms in 
random potentials have been reported For trapped 
cold atoms, a lattice potential can be implemented by far 
detuned counter-propagating laser beams, the AC driving 
can be obtained using a periodic phase-shift between the 
beams |32| and the random potential can be obtained with 
a superimposed random speckle pattern [23 ■ We believe 
our results could lead to very interesting effects in these 
systems. 
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